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We study the production of gravitational waves during oscillations of the inflaton around the 
minimum of a cuspy potential after inflation. We find that a cusp in the potential can trigger copious 
oscillon formation, which sources a characteristic energy spectrum of gravitational waves with double 
peaks. The discovery of such a double-peak spectrum could test the underlying inflationary physics. 


Introduction. Gravitational waves play an important 
role in the context of inflationary cosmology. A stochas- 
tic background of gravitational waves, produced during 
inflation and subsequent preheating/reheating after in- 
flation, carries useful information about the inflationary 
dynamics and inflaton decay (see |1] for a recent review). 
Detecting such a stochastic background of gravitational 
waves, whether directly or indirectly, can provide us with 
a unique opportunity to test theories of inflation. 


During inflation, quantum fluctuations of the tensor 
modes of the spacetime metric were stretched by the ac- 
celerated expansion of the Universe, and were then nearly 
frozen on super-Hubble scales. In the standard single- 
field slow-roll inflationary scenario, the amplitude of the 
power spectrum of the resulting gravitational waves de- 
pends on the energy scale of inflation. Since these grav- 
itational waves can result in B-mode polarization of the 
cosmic microwave background (CMB) anisotropies, their 
spectrum is in principle measurable by CMB polariza- 
tion experiments. Although primordial B-modes have 
not been detected yet, the amplitude of the tensor spec- 
trum (quantified by the tensor-to-scalar ratio r) has be- 
come a strong discriminator of models. Current CMB 
data alone already put_an upper bound on r < 0.09 at 
95% confidence level fat, and when combined with the 
constraints on the scalar spectral index, have been effec- 
tive in discriminating inflationary models. For example, 
the cubic and quartic potentials are strongly disfavored, 
and the quadratic potential is moderately disfavored by 
the Planck 2015 data EN while axion monodromy infla- 
tion with a linear potential Zi or fractional powers E are 
compatible with the current Planck results. Further ad- 
vances in axion monodromy inflation have suggested po- 
tentials with even more possible powers |6, (7. Moreover, 
it has recently been shown that stringy effects can lower 
the power of a quadratic axion monodromy potential to 
less than linear Thus, axion monodromy inflation 
represents an interesting class of large field inflationary 
models that are compatible with data. 


Besides vacuum fluctuations during inflation, another 
source of gravitational waves is parametric resonance 
during preheating (9|. During preheating after inflation. 


the Fourier modes of a scalar matter field x coupled 
to the inflaton grow exponentially by parametric reso- 
nance, driven by the oscillating inflaton. The resonant 
modes are quickly pumped up to a large amplitude. Such 
highly pumped modes correspond to large, time depen- 
dent density inhomogeneities in position space, ensuring 
that the matter distribution has a non-trivial quadrupole 
moment, which can source a significant spectrum of grav- 
itational waves [10]. The present peak frequency of such 
gravitational waves is proportional to the energy scale of 
inflation {14}, while the present amplitude of gravitational 
waves is independent of the energy scale of inflation ug. 
In hybrid inflation with the energy scale ranging from 
the GUT scales down to the TeV scale, the stochastic 
background produced during preheating is expected to 
be directly detected by future gravitational wave detec- 
tors ; 

In this letter, we shall investigate the production of 
gravitational waves during oscillations of the inflaton af- 
ter inflation with a cuspy potential 


V (e) = AMS "lel? (1) 


with p = 1, 2/3, 2/5, and Mp1 = (87G)~!/? is the reduced 
Planck mass. Our study is motivated by the potentials 
that appear in axion monodromy inflation [4 Bi. though 
we hasten to add that the powers of these potentials are 
expected only at large field values, due to the coupling of 
the inflaton to high scale physics. At the end of inflation, 
i.e., for small ¢, these potentials for axion monodromy be- 
come quadratic. Nonetheless, cuspy potentials can arise 
in other inflationary contexts, e.g., through non-standard 
kinetic terms or as a result of integrating out the dynam- 
ics of other fields that couple to the inflaton. Thus, we use 
these specific potentials as benchmarks to illustrate the 
novel gravitational wave signatures that can arise when 
the potential has a cuspy behavior at the end of inflation. 
Assuming the potential in eq. applies to both the in- 
flationary era and at the end of inflation, the value of À in 
this simple class of models can be fixed by the estimated 
amplitude of scalar perturbations from the CMB data. 
For powers of p = 1, 2/3, 2/5, A + 3,4,5 x 10719, the pre- 
dicted scalar spectral index n, ~ 0.970, 0.973, 0.976, and 
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the predicted tensor-to-scalar ratio r ~ 0.08, 0.05, 0.03, 
respectively, assuming the number of e-folds N = 50. 
These predictions are in agreement with the recent CMB 
data. In the reheating scenario, the inflaton ¢ oscillates 
near the minimum of its potential after inflation and de- 
cays into elementary particles. However, due to the cusp 
of the potential, the oscillating behavior of the inflaton is 
very different from that of smooth potentials like ¢? and 
#*. It has recently been shown that an efficient paramet- 
ric resonance can occur during preheating for an inflaton 
potential of the form of eq. (i) with 0 < p < 2, if the 
inflaton is coupled to a scalar matter field y via an in- 
teraction term $7? (14). However, the production of 
gravitational waves has not been studied to our knowl- 
edge. In this letter, we are interested in gravitational 
wave production during oscillations of the inflaton after 
inflation with cuspy potentials. We find that the non- 
smooth oscillations can trigger amplification of fluctua- 
tions of the inflaton itself at the moment when ¢(t) = 0, 
so that oscillons copiously form after inflation. As in the 
models with a symmetric smooth potential and an 
asymmetric smooth potential lid), the oscillon formation 
in the models with cuspy potentials sources a stochastic 
background of gravitational waves, on which the charac- 
teristic size of the oscillons is imprinted. Interestingly, 
these cuspy potentials result in a characteristic energy 
spectrum of gravitational waves with double peaks, which 
can be distinguished from smooth potentials by probing 
the shape of the energy spectrum of gravitational waves. 
Model. Gravitational waves are described by the 
transverse-traceless gauge-invariant tensor perturbation, 
hij, in a Friedman-Robertson-Walker (FRW) metric, 


ds? = —dt? + a° (t)(ôij + hij)dx'da’, (2) 


where t is the cosmic time and a(t) is the scale factor. To 
first order in hij, the perturbed Einstein equation reads 


2 
—— JTT 
Ma’ ay? (3) 


hij + 3Hhij — Vhs = 
where mt is the transverse-traceless projection of the 
anisotropic stress tensor T;j. In our model we assume 
that the inflaton is weakly coupled to other fields during 
preheating. Gravitational waves are sourced by the infla- 
ton fluctuations 64, i.e., IT = (0;560;5¢)"". Later we 
shall discuss the effects of interactions with other fields 
on the production of gravitational waves. The energy 
density of gravitational waves is given by 


= Mag, i 4 
PGW = FT ij ) ( ) 


where (...) denotes a spatial average over the volume. It 
is commonly parameterized by the dimensionless density 
parameter per logarithmic frequency interval, QNqwh? = 
h?dpaw /dink/pc, where pe is the critical density of the 


Universe and h is the reduced Hubble parameter. 

To understand oscillon formation during oscillations 
of the inflaton after inflation with a cuspy potential, we 
now investigate the evolution behavior of inflaton fluc- 
tuations. To first order, the equation of motion for the 
Fourier modes of fluctuations of the inflaton ¢ is 


2 


Sbr + 3Hõbr + = + aC) Ör =0. (5) 


To solve this equation, we neglect the expansion of the 
Universe for the moment, thus the friction term drops out 
of the equation of motion. Moreover, since we are inter- 
ested in large-scale modes, the k? term can be dropped. 
The equation becomes 


don + V"($)5bx = 0. (6) 


For illustrative purposes, in what follows let us consider 
the linear potential. We shall later numerically show that 
there exist similar evolutions of 66, for p = 2/3 and 
p = 2/5. Since the inflaton potential V(¢) = AM3|¢| 
has a cusp at ¢ = 0, its derivative with respect to ¢ is 
a step function and its second order derivative is a delta 
function V”(ġ) = 2\M3,6(¢). We focus on the evolu- 
tion behavior of the y modes near the point ¢(t) = 0. 
It is convenient to define t such that d(t = 0) = 0. 
The solution to the equation of motion for the infla- 
ton is O(t) = |dm|t + AM3t?/2 when t < 0 and ¢(t) = 
ideale AM3t?/2 when t > 0, where Èm is the max- 
imum value of ¢ at t = 0. Since ¢(t) ~ |dmlt is a 
good approximation in a small vicinity of ¢ = 0, we find 
ddg(t = 0*) = ddg(t = 07) = —2AM3 dor (t = 0)/lml, 
which implies that 6¢, jumps suddenly when ¢ crosses 
the cusp of the potential. Such periodic jumps of Sk 
result in periodic, rapid increases of the energy density 
pr. We show the time evolution of the energy density 
for a linear potential (orange) in an expanding Universe 
in Fig. [I] In the cases of p = 2/3 and p = 2/5, |V" (¢)| 
becomes infinite when |¢| tends to zero. To avoid this 
singularity, the potential with a cut-off |¢| > O(107%) is 
adopted in our numerical calculations. We also show the 
time evolution of the energy density for the ¢?/° poten- 
tial (blue) and ¢?/° potential (green) in Fig. [I] Similar 
to the linear potential, the energy density suddenly in- 
creases near the points at which $(t) = 0. However, the 
sudden increase is followed by a sudden decrease near 
g(t) = 0 after several oscillations of the field ¢. The 
energy density always increases after each oscillation of 
the field ¢. We have checked that the increase in energy 
density is independent of the choice of the cut-off. As a 
result, oscillon formation occurs when the inflaton oscil- 
lates near the minimum of its potential, which sources a 
stochastic background of gravitational waves. 
Simulation Results. Using a modified version of LAT- 
TICEEASY aA, we simulate the production of gravita- 
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FIG. 1: Evolutions of the energy density pẹ for cuspy poten- 
tials with p = 1 (orange), p = 2/3 (blue) and p = 2/5 (green), 
respectively. 


tional waves during preheating in the models (i) with 
cuspy potentials. The publicly available C++ package 
LATTICEEASY has been developed for more generally 
calculating the evolution of interacting scalar fields in an 
expanding Universe using a staggered leapfrog algorithm, 
which allows us to simultaneously solve the Friedmann 
equations and the inflaton evolution equation in an FRW 
Universe. To calculate the energy spectrum Qqwh? of 
gravitational waves, the gravitational wave equation (B) 
needs be solved in an expanding Universe. The spectral 
method can directly solve for the evolution of hij sourced 
by the transverse-traceless part of the anisotropic stress 
tensor T;j in Fourier space iş. Actually, one can first 
evolve the tensor perturbations in configuration space 
sourced by the anisotropic stress tensor T;j instead of its 
transverse-traceless part and then apply the transverse- 
traceless_projector to the real physical hj; in Fourier 
space ag. Another method is based on the Green’s 
functions in Fourier space to calculate numerically the 
energy spectrum of gravitational waves generated well in- 
side the horizon . In our lattice simulations we adopt 
the configuration-space method for solving the following 
evolution equation of the tensor perturbations 


2 


: ; 1 9 
tig +3Hüúi; — zy Uij = Ma? 


Tr (7) 


Therefore, the transverse-traceless tensor perturbations 
can be written as h;;(t,k) = Mij.tm(k)uim (t, k), where 
Nij im (k) is the transverse-traceless projection operator 
and Ujm(t,k) is the Fourier transform of the solution to 
the equation (Z). The energy density of gravitational 


waves can be expressed in terms of u;j as 


paw = J5 [ata Big (1) (tK). (8) 


We perform three-dimensional lattice simulations with 
256° points in a box with periodic boundary conditions 


assuming À = 1.26 x 1071? in the linear potential model. 
We set the initial values of the inflaton, its derivative 
and the scale factor as $; = 0.75Mpi, ĝi = 6.8 x 1074M 
and a; = 1. The inflaton fluctuations and its derivative 
are initialized by quantum vacuum fluctuations, while 
the tensor fluctuation and its derivative are initialized as 
zero. We stop the simulation when the energy spectrum 
of gravitational waves does not grow significantly. The 
energy spectrum and its frequency at the end of simula- 
tions are converted to the present values. Fig. [2] shows 
the time evolution of the energy density as a function of 
position on a two-dimensional slice through the simula- 
tion from a(t) = 6.13 (top-left), to 6.81 (top-right), to 
7.43 (bottom-left) and 13.4 (bottom-right) in the linear 
potential model. We can see that at the beginning, os- 
cillons copiously form and then decay during oscillations 
of the inflaton. In this model, the rapid growth of oscil- 
lons results in the production of gravitational waves with 
Qawh? ~ 2 x 107° today. Our lattice simulation results 
show that there appear two peaks in the energy spec- 
trum of gravitational waves, a feature very distinct from 
that of other models. Therefore, our model can be distin- 
guished from the production of gravitational waves dur- 
ing preheating by future gravitational wave detectors. A 
phenomenological study of gravitational waves produced 
from preheating with a time dependent resonance param- 
eter q(t) was recently undertaken a. For some choices 
of q(t), one also finds a gravitational wave spectrum with 
multiple peaks due to non-linear effects. Our model pro- 
vides an explicit realization of this phenomenon. 

The evolution of the spectrum goes through two differ- 
ent stages, the linear growth stage and nonlinear growth 
stage. In the first stage, as shown in Fig. [3| the small-k 
modes of the field ¢ exponentially grow due to the cusp 
of the potential until the turning point a(t) = 7.12. The 
linear growth is more efficient than those driven by the 
symmetric potential and asymmetric potential ig. 
This leads to the left peak in the energy spectrum of 
gravitational waves, which is characteristic of the cuspy 
potential. In the second stage, from Fig. B] we see that 
the small-k modes begin to drop and the large-k modes 
continue to grow. It implies that the energy flows from 
the small-k modes to the large-k modes, as discussed in 
detail in Eg. This leads to the right peak in the energy 
spectrum of gravitational waves. 

Moreover, from the lattice simulations of preheating in 
the models (I) with p = 2/3 and p = 2/5, we find that 
the energy spectra of gravitational waves peak at around 
Qawh? ~ 1.2 x 107° and Qawh? ~ 4 x 10710, respec- 
tively, which are lower than the one in the linear potential 
model. Although the energy density pp increases more 
than in the p = 1 case when |¢| tends to zero, as shown 
in Fig. [L| a sudden decease follows, which may suppress 
the production of gravitational waves. 

In our analysis, we have neglected the interactions be- 
tween the inflaton ¢ and other matter fields. If the in- 
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FIG. 2: Energy density p on a two-dimensional slice through the simulation, when a(t) = 6.13 (top-left), 6.81 (top-right), 7.43 
(bottom-left) and 13.4 (bottom-right), in the linear potential model. 
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FIG. 3: Evolution of the energy density of the field ọ for 
the linear potential. The yellow line corresponds to a turning 
point a(t) = 7.12. 


flaton is coupled to a matter field x, broad parametric 
resonance leads effectively to a fast growth of the fluc- 
tuations of x 14]. However, our numerical simulations 
confirm that the growth of the inflaton fluctuations them- 
selves triggered by the cusp in its potential is more ef- 
fective than that of the field y by parametric resonance. 
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FIG. 4: Energy spectra of gravitational waves today, pre- 
dicted by the linear potential with A = 7.5 x 107*8 (blue) 
and À = 9.72 x 107*° (green). The black curve is the ex- 
pected sensitivity curve of the fifth observing run (O5) of the 
aLIGO-Virgo detector network. 


Therefore, gravitational waves are sourced mainly by the 
inflaton fluctuations, even if a parametric resonance for 
the field x occurs in this model. 

Observational Implications. As found in (12), Ep the 
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peak amplitude of the energy spectrum of gravitational 
waves needs not depend on the energy scale of inflation, 
while the peak frequency scales inversely with the energy 
scale of inflation. In the single-field slow-roll inflation- 
ary scenario, if A ~ 3 x 1071 is fixed by the amplitude 
of the primordial power spectrum of curvature perturba- 
tions A, = 2.2 x 107°? [3], the peak frequency of gravita- 
tional waves today is fixed to be f ~ 10°Hz, many orders 
of magnitude beyond the frequencies that can be reached 
by current gravitational-wave detection experiments. 

If the model parameter is not fixed by the amplitude 
of the primordial power spectrum of curvature pertur- 
bations, the sensitivity of advanced LIGO (aLIGO) is 
expected to be significantly improved, which allows us to 
possibly observe gravitational waves produced during os- 
cillations of inflaton after inflation. For example, in the 
hybrid inflationary scenario (21), since ¢ is not necessarily 
the inflaton itself, A becomes essentially a free parame- 
ter. In this case we have plotted in Fig. 4]the present-day 
energy spectra of gravitational waves produced during 
oscillon formation in the linear potential model with 
A = 7.5 x 10743 (blue) and A = 9.72 x 1074% (green), the 
peaks of which lie above the expected sensitivity curve of 
the fifth observing run (O5) of the aLIGO-Virgo detec- 
tor network BRA. From Fig. [4] we can see that there are 
two peaks in the energy spectrum of gravitational waves, 
which differ from other spectra of gravitational waves 
produced during preheating. A detection of the second 
peak may require corroboration from other gravitational 
wave detectors such as the Big Bang Observatory (BBO). 

To summarize, we have studied the production of grav- 
itational waves during oscillon formation after inflation 
with cuspy potentials. At the end of inflation, oscil- 
lon formation can be triggered by the particular oscilla- 
tions of the inflaton around the minimum of its potential, 
which sources a characteristic double-peak spectrum of 
gravitational waves. The discovery of such a background 
would open a new observational window into inflationary 
physics. 
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